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ABSTRACT 





In structural redesign, two structural states are involved; the baseline (known) State SI with 
unacceptable performance, and the objective (unknown) State S2 with given performance 
specifications. The difference between the two states in performance and design variables may be 
as high as 100% or more depending on the scale of the structure. A Perturbation Approach to 
Redesign (PAR) is presented to relate any two structural states SI and S2 that are modeled by the 
same finite element model and represented by different values of the design variables. General 
perturbation equations are derived expressing implicitly the natural frequencies, dynamic modes, 
static deflections, static stresses, Euler buckling loads and buckling modes of the objective State 
S2 in terms of its performance specifications, and State SI data and FEA results. LargE 
Admissible Perturbation (LEAP) algorithms are implemented in code RESTRUCT to define the 
objective State S2 incrementally without trial and error by postprocessing FEA results of State SI 
with no additional FEAs. Systematic numerical applications in redesign of a 10-element 48-d.o.f. 
beam, a 104-element 192-d.o.f. offshore tower, a 64-element, 216-d.o.f. plate, and a 144 element 
896-d.o.f. cylindrical shell show the accuracy, efficiency, and potential of PAR to find an 
objective state that may differ 100% or more from the baseline design. 

I. INTRODUCTION 


Several problems in analysis, design, and modification of a structure or a structural design 
can be stated as redesign problems. Those are two-state problems involving the baseline State SI 
and the objective State S2. SI is known and has been modeled and analyzed by FEM. In the 
event that the performance of State SI is unacceptable, the objective State S2 must be defined to 
satisfy performance specifications. The Perturbation Approach to Redesign (PAR) developed in 
this work can relate any two structural states that can be modeled and analyzed by the same FE 
model. PAR has the potential to perform redesign in the sense of resizing, reshaping, and 
reconfiguration to satisfy any performance requirements that can be predicted by FEA including 
modal dynamics, static deflections and stresses, and global buckling. LEAP algorithms 
implemented in code RESTRUCT (REdesign of STRUCTures) [3] presently can handle resizing 
for natural frequencies, mode shapes and static deflections. 

Figure 1 shows several two-state problems that appear in the analysis-design-redesign 
process following a basic FE analysis. In analysis, the following two-state problems are 
encountered: (PI) Model correlation [28], (P2) Derivation of global failure equations [1, 14], (P3) 
Failure point identification [14], (P4) Redundancy [14], (P5) Reliability, [4], (P6) Non- 
Destructive-Testing [24]. In design, the following two-state problems are encountered: (P7) 
Redesign for target performance [1, 2, 11, 12, 24, 26, 27], (P8) Redesign for target redundancy, 
(P9) Redesign for target reliability. 

LEAP theory was developed during the past seven years from the linear perturbation 
techniques introduced by Stetson in 1975 [26, 27] and modified by Sandstrom et al [24], They 
redesigned a structure for both natural frequency and mode shape objectives but allowed only 
small differences between the baseline and objective states. In that respect, linear perturbation 
methods are equivalent to design sensitivity methods. Nonlinear perturbation methods [11, 12] 
allow for large differences between the two states. The objective state is found by postprocessing 
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data of the baseline structure only, using an incremental prediction correction scheme [1], 
Presently, research efforts are directed towards two goals. The first one is to redesign larger scale 
structures as far away from the baseline structure as possible before a second FEA is needed. 
Large admissible perturbations [1] updating only cognate modes [2] in an incremental process are 
used towards that end. Substructuring is also investigated for that purpose, as well as for 
reshaping and reconfiguration. The second goal is to implement more and different objectives and 
derive the corresponding general perturbation equations. LEAP algorithms are under development 
for static stress, global buckling load, and buckling mode objectives. 

The problem of redesign by large admissible perturbations is analyzed in Chapter II. Several 
two-state problems mentioned above are stated as redesign problems in Section ILL The 
Perturbation Approach to Redesign (PAR) is presented in Section II. 2 and LEAP theory for 
development of solution algorithms is summarized in Section II. 3. Many numerical applications 
using four different structures are presented in Chapter III to assess the present status of code 
RESTRUCT, and the potential and limitations of PAR. 

II. REDESIGN BY LARGE ADMISSIBLE PERTURBATIONS 

A simple modeling-analysis-design-redesign process for structures using FEM is shown in 
Figure 1. Rectangular blocks indicate two-state problems which can be formulated as redesign 
problems using PAR and solved efficiently by a LEAP algorithm. Shaded blocks indicate 
problems already solved in some form by code RESTRUCT. Some of those problems are 
discussed below. 

II.l. Redesign and Other Two-State Problems 

The classical structural redesign problem appears in Figure 1 after analyzing either the original 
or the correlated FE model. Undesirable response - such as a natural frequency in the range of 
wave excitation, a dynamic mode with high amplitudes near the free surface where wave and 
current loads are maximum, or high stresses and deflections - makes redesign mandatory. The 
performance specifications of the objective design are desirable values of those response 
particulars. 

After placing a structure in service, tests are performed to measure its performance and 
compare it to FEM predictions. In the modeling process, simplifying assumptions, uncertainty, 
and ignorance result in discrepancies between measurements and predictions particularly for 
marine structures which have large manufacturing tolerances. The process of finding a FE model 
of a physical structure that will correctly predict measured structural response is called model 
correlation. The initial FE model is the known State SI. The objective State S2 represents the 
unknown correlated FE model. The Perturbation Approach to Redesign presented in the following 
section preserves element connectivity and changes geometric properties so that the correlated 
model represents a real structure [2]. That is, PAR does not change simply numbers in the mass 
and stiffness matrices. PAR can also solve the problem of model correlation for geometry 
dependent hydrodynamic load [28]. 

The problem of failure point identification can also be formulated by PAR and solved by a 
LEAP algorithm. SI represents the initial structural state and S2 the unknown failure point 
(design point in reliability terminology) on a limit surface [10, 20]. The advantage of PAR is that 
it can provide an implicit expression for a global failure criterion by relating State S2 to SI. 

Related is the problem of reserve and residual redundancy. In the literature, several different 
aspects of redundancy are presented as definitions depending on the type of structure and analysis 
performed [5, 6, 21]. PAR remedies this lack of invariant and consistent redundancy definition by 
introducing a redundancy injective mapping [14] defining the difference between the initial intact 
or damaged structure and the design point 

Finally, a new methodology for reliability analysis and design of large scale structures is 
under development based on PAR [4], The Perturbation Approach to Reliability provides an 
alternative to the systems approach [5, 21, 29] and the stochastic FEM [19, 30] which are the two 
most popular methods in structural reliability. PAR makes possible the introduction of advanced 
structural analysis in the reliability computations without simplifying the structure. PAR also 
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T allows randomness in geometry, material, and load. There is no limitation to the number of 
random variables used and the random load need not be applied incrementally until structural 
failure. The reliability analysis problem is a two-state problem where S 1 is the initial structure and 
S2 the design point. 

PAR can also address the very difficult problems of target redundancy and target reliability 
design. SI is the initial structural design of inadequate redundancy or reliability and S2 is the 
objective structure of specified redundancy or reliability [4], These are difficult design problems 
because redundancy and reliability are not computed by FEM. PAR can solve these problems 
because of the introduction of an injective mapping relating S 1 to S2. 


II. 2. Perturbation Approach to Redesign (PAR) 

The PAR methodology has been developed to solve the above two-state problems. It has five 
major steps: Step 1 : A Structure (SI) is modeled and analyzed by a general purpose FE code; 
MSC/NASTRAN is used in our work. So far, four types of analysis have been considered in 
PAR and the governing equations are listed below. For modal dynamics the free vibration 
equations for S 1 are 

([k] - (Oj [m]){v}j = {0} for j=l,2,...,n , (1) 

where then eigenvalues coj , j=l,2,...,n satisfy equation det ([k] — coj [m]j = 0 . Inequation 

(1), damping may be included only in Rayleigh's form and added mass is included in [m] . For 
the static deflections and stresses of SI we have 


[k] {u} = {f} (2) 

and {cr} = [S] [k] 1 {f} , (3) [S] = [G] [D] [N] , (4) 

where [G], [D], [N] are the stress-strain, strain-displacement, and shape function matrices. Th 
governing equation for global buckling in finite elements is 

(M + MWbMO} . < 5 > 

where [Iq,] and [k a ] are the small displacement and initial geometric stiffness matrices. 

Step 2 : The following perturbation relations are introduced relating State S2 to S 1 : 

[k'] = [k] + [Ak] , (6) [m'] = [m]+[Am] , (7) 

rco'2^1 = r~a> 2 -J + rA(o> 2 )-J , (8) [<t>’l = [<t>] + [A<>] , (9) 


where unprimed and primed symbols refer to the baseline (SI) and the objective State S2, 
respectively, and prefix A indicates difference between counterpart quantities of states S 1 and S2. 


L 


W = [{¥},. {\|r} 2 , ...,{y} n ] , is the matrix of eigenvectors of SI and ['co 2 ~J is the diagonal 
matrix of the corresponding eigenvalues. Perturbation relations pertaining to equations (2) and (3) 
are 

{u*} = (u) + {Au} , (10) 

{a'} = {cr} + (Aa) , (12) 

For the global buckling eigenvalue problem we have 
[ki] = [ko] + [Ak 0 ] , (14) 


{f '} = { f } + { Af } , 

(id 

[S'] = [S] + [AS] . 

(13) 

[k a ] = [k c ] + [Ak a ] , 

(15) 


J 
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r [P c 'r] = [Pcr] + [AP cr ] , (16) K] = K] + [Ad> b ] • (17) 1 

Further, in Step 2, desirable values of some response particulars of S2 such as natural 
frequencies, static deflections and mode shapes are specified. An incomplete set of mode Shapes 
may be used and only some degrees of freedom may be defined in each mode. 

Step 3 : The differences in structural properties between SI and S2 are expressed in terms of the 
fractional changes a e ,e=l,2,...,p of p properties of elements or groups of elements as: 

[Ak] = X[ Ak e] = L[ k e] a e . (18) [H= X[ Ara e ]= IKK .(19) 

e=l e=l e=l e=l 

[AS]=XKS e ] = X[ S eK • (20) 

e=l e=l 

Several a e s may refer to the same element but different properties such as bending, torsion, and 

stretching. The unknowns in the process of defining S2 from its specifications and S 1 are the 

fractional changes a e . When the a e s are defined it is ensured that element connectivity in the 
FE model is preserved and S2 represents a real structure. 

Step 4 : The differences in structural response between states SI and S2 are expressed implicitly in 
terms of the o^s by the general perturbation equations. For modal dynamics we have 


KeKv'L - K]{v'}i 

e=l 

ZMj[ke]{v'} i a e =-{v'}[[k]{v'} i - (22) 

e=l 

£{\ j /'}J[m e ]{v , } i a e =-{v , }J[m]{\j/’} i , (23) 

e — 1 

for i=“I,2,„.,n, j = i+1, i+2,...,n [1,2] . Equation (21) represents the n diagonal terms of 
the energy balance equation = 0 for S2, that is, the Rayleigh quotients for 

CDp . Equations (22) and (23) represent the orthogonality conditions of modes with 

respect to [k'] and [in'] . Theoretically, orthogonality of modes with respect to one of [k’] or 
[m'] implies orthogonality with respect to the other. Numerically, however, both conditions must 

be forced if [\|/'}j, j=l,2,...,n, are to represent modes of a real structure. 

The general perturbation equations for static deflections are derived from the counterpart of 
equation (2) for structure S2 based on the modal dynamic expansion of {u'} in terms of the 


)a e = (D[ 2 {\(f'}j r [m]{\i/'} i -{v'}j r [k]{\|/'} i , (21) 



unknown modes {y'jjJ = l,2,...,n . Thus, inversion of matrix [k'] is avoided. Linearizing 
only the explicit dependence on the a e s, we have [1,15] 


_ V f ^im A m 
B r 


u i = L 


4*im 


me 


m-1 ' 
n 


'm 


a 


e » 


where 


-I 

' e=lVm=l 

={ V }^[k e ]{ V } . 
j=l 

The general perturbation equations for static stresses are derived in a similar manner [14] 


(24) 


L 


{Aa} = -{a} + 


[S]+ SfSefo 

e=l 


n x/ 


a' a P f n , / A 
V 9 im A m _ Y 

Zj k jL Li ~ n 2 ^rn e 
J |_m=l m e=lVm=l 


a P 


(25) 
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~For global buckling, the general perturbation equations are derived using the same method as in j 
the case of the modal dynamics eigenvalue problem [14] 


i M,}’' (K ] - nK, ]){vb}i« e = Mi 1 ' (ivK ] - M m, • 

e=l 

(26) 

e=l 

e=l 

(27) 

(28) 

for i = l,2,...,n , j=i + l,i + 2,...,n , where [k c ] = [k 0 ] - [l^p] , k oF includes the body force, 
and [kd = - Pi[k CT0 ] - [k cF ] . 

Step 5: In this final step, the problem of finding State S2 based on its specifications and results of 
FEA for SI is formulated and solved for the p unknown a e s using the LEAP algorithm 
presented in the next section. The problem formulation is as follows: 

Minimize liaise 91 p , 

(29) 


subject to n^ natural frequency objectives co[ 2 , i - 1,2,. ...n^ ; n^ normal mode objectives 

fa , number of (k, i) = n^ ; n u static deflection objectives u- , i=l,2 n u ; n CT static stress 

objectives o- , i=l,2,...,n a ; n b global buckling eigenvalues P( , i=l,2,...,n b ; n^ b buckling 
mode objectives , number of (k,i) = n^ b ; 2p lower and upper bounds on the redesign 

variables ote.-K a; <a e <a^, e=l, 2 p; n a admissibility constraints extracted from 

equations (22) and (23), where n a = 2^(n r -i) = n^^n,. -^-n^j ; and n ab 

i=l 

admissibility constraints extracted from equations (27) and (28), where 

n ab = 2 S( n r “ 0 = n b[( 2n r -1)- "b] • Ml of *** above objectives are substituted in 

the appropriate general perturbation equations (21)-(28). The remaining unused general 
perturbation equations may be used to predict the unspecified performance particulars of the 
objective State S2. Accuracy of those predictions, however, it not as high as those of the redesign 
objectives. All the constraints of the above problem may result in an empty, non-empty, or 
countable feasible domain. In the first case, the redesign objectives cannot be achieved for the 
selected set of redesign variables, in which case a minimum error solution in satisfaction or the 
redesign objectives is achieved by a generalized inverse algorithm [1, 2, 11, 15]. In the second 
case, an optimum solution is achieved using an optimality criterion (29). 

II. 3. T .arpF. Admissible Perturbation (LEAP) Algorithm 

The redesign problem formulated by PAR in Section n.2 can be solved by a LEAP algorithm. 
Many LEAP algorithms have been developed to solve a variety of two-state problems [1, 2, 14, 
15 28] and have been documented in detail. Suffice to present here the basic steps and difficulties 
of ’the solution algorithm. The LEAP algorithm developed to solve the redesign problem is 
outlined in Figure 2. It starts from the baseline structure (SI) and reaches incrementally thej 
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f objective S2 by prediction and correction. In the prediction phase of the algorithm, the small ~j 
I perturbation method [24, 26, 27] is used. The modal dynamics general perturbation equations are I 
linearized. For that purpose, increments are limited to 7% differences between S2 specifications 
and the corresponding SI properties. Predictions are small but inadmissible because admissibility 
conditions (22) and (23) are linearized. In the correction phase , perturbations are corrected by 
satisfying the nonlinear general perturbation equations and are forced back into the admissible 
space by satisfying the nonlinear admissibility conditions. The total CPU time for redesign may 
be reduced by a factor occasionally as high as 4 when in the first increment the space of cognate 
modes is identified and thereafter all computations are performed in that space. Such is the case 
for torsional redesign [2] of the offshore tower in Figure 5. Torsional modes (3, 18, 19) 
constitute one cognate subspace with very weak interaction with other modal subspaces such as 
those for bending and stretching. 

In each increment, in both phases the resulting problem may be underdetermined or 
overdetermined depending on the relation between the number p of redesign variables a e , the 
number of equality constraints (S2 specifications) n = n ffl + n^ + n u + n CT + nj, + n^ + n a + n a ^ , 
and the 2p bounds on the a e s. When the problem is overdetermined, a minimum error solution 
in satisfaction of the S2 specifications is produced by a generalized inverse algorithm. When the 
problem is undetermined, it is solved by optimization using the minimum change criterion in 
equation (29). To achieve this global objective, at each increment the following objective is 
minimized 

p r /-i 

min £ (l+/a e )Il(l + q a e )-l • (30) 

e=l 4=1 

The problem is solved by quadratic programming [8] or sequential quadratic programming [7] 
depending on whether the expression for [Ak] is linear as in equation (18), or nonlinear as in the 
case of plate and shell redesign. In those cases, the plate or shell thickness is selected as redesign 
variable resulting in a cubic expression for [Ak] in terms of the c^s . [AS] is always a nonlinear 
expression of the a e s because [S e ] depends on the distance of the point where the stress is 
computed from the neutral axis. The LEAP algorithm is implemented in code RESTRUCT 
(REdesign of STRUCTures) [3]. It is 27,000 FORTRAN 77 commands and may serve as a 
postprocessor to any special or general purpose FE code. We presently use it to postprocess 
MSC/NASTRAN.V64 data on the secondary (UB) main frame computer (IBM-3090) of the 
University of Michigan. 

The LEAP algorithm outlined above finds the optimum objective structure S2 without trial and 
error and with no additional FEAs. The [k] matrix inversion required in static deflection and stress 
redesign is avoided by using modal expansions as shown in equations (24) and (25). Thus, an 
accurate modal basis is mandatory even as S2 moves far away from SI. LEAP algorithms can 
surmount the following three difficulties as well. All general perturbation equations (which 
become equality constraints in the optimization problem) are strongly nonlinear implicit 
expressions of the redesign variables a e . The static force vector [f } may depend on the 
structure's geometry (e.g. hydrodynamic loads) and consequently change in the redesign process. 
Finally, the set of specifications provided for S2 are usually incomplete and only some d.o.f.s of 
specified modes are defined. 

III. NUMERICAL APPLICATIONS 

A total of 42 numerical applications are presented in this section on optimal redesign of four 
different structures [22, 9, 31]. Results are summarized in Tables 1, 4, 5, 6 and show the 
accuracy of code RESTRUCT for applications with number of redesign variables ranging from 8 
to 21; natural frequency and mode shape redesign objectives changing by a factor ranging from 
0.3 to 2.0; degrees of freedom ranging from 48 to 896. For each redesign objective, Tables 1, 4, 
5, 6 show the objective value, the value actually achieved as computed by reanalysis with 
MSC/NASTRAN and the corresponding relative error. CPU time and numbers of extracted 
■ modes n r , admissibility conditions n a , and redesign variables are also shown. The values of the . 
[_ redesign variables of the optimum solution are not shown. The optimal solution appears in thej 
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form of optimal Eucledian norm of the a e s in Tables 5 and 6; and in the form of the Hasover-] 
Lind reliability index [10] in Tables 1 and 4. 

10-element 48-d.o.f. beam : The clamped-hinged beam in Figure 4 is subjected to a uniform load 
in the y direction and a concentrated load applied at node 7 in the z direction. CO! = 183.092 
rad/sec, the horizontal and vertical deflections at node 7 are v-j = 12.151 mm and w 7 = 17.7o3 
mm as computed by MSC/NASTRAN. Redesign variables and structural groups are shown in 
Table 2. The accuracy of the redesign process is shown in Table 1 for one, two or three 
simultaneous redesign objectives. The problem of reliability analysis is studied assuming 
randomness in geometric properties, A (area), I y , I z (moments of inertia), and material properties 
E (Young's modulus) and p (density). The fractional changes a e are assumed to be 
independent normal random variables of zero mean. Standard deviations are selected as CJ a£ j = 
0.40 for bending rigidities EI y and EI Z , and c apA = 0.30 for mass per unit length pA . In 

order to compute the probability of failure to first (FORM) or second (SORM) order [20], 
computation of individual and joint design points and the corresponding Hasover-Lind reliability 
index (3 is required as shown in Figure 3. Computation of fi is achieved by transforming the 
OgS to independent standard normal random variables through the Rosenblatt transformation [13]. 
These numerical applications as well as those following on the offshore tower show that large 
admissible perturbation methods can introduce sophisticated structural analysis in reliability 
without simplifying the structural model and without repeated FEAs [4], 

104-element 192-d.o.f. offshore tower : The offshore tower shown in Figure 5 is 69.95 m high 
and operates in 45.72 m water depth. The tower at the base is square with a 38.10m side and 
tapers linearly to 22.86 m at the deck. The FE model of the tower is composed of 104 circular 
tubular beam elements and has 192 dofs. Loading on the tower is due to: (i) 240 tonnes deck load 
which is applied to the structure as uniformly distributed load at the deck nodal points, (ii) Wave 
hydrodynamic forces calculated for a design wave of 182.88 m length and 6.10 m height using 
Morison's equation. The wave propagates in the x-direction. (iii) Wind generated water current in 
the x-direction with linear velocity profile of 1.03 m/sec at the mean free surface waterline and zero 
at the sea bed. coj = co 2 = 4.695 rad/sec for the first bending modes in the XZ and the YZ 
planes. 003 = 5.353 rad/sec for the first torsional mode with respect to axis Z . Redesign 

variables and structural groups are shown in Table 3. 

Failure states are defined by deterioration factors in the first and third eigenvalues of 1.54 and 
2.00. Geometric and material properties are random. The fractional changes a e s , shown in 
Table 3, are assumed to be independent normal random variables with zero mean. Standard 
deviations are selected as o a£1 = 0.40 for bending rigidity El and a ttpA = 0.30 for mass per 

unit length pA . Design points are again computed by postprocessing FE analysis results for the 
baseline design only. It should be noted that both in Tables 1 and 4 the computed p are very high 
because the external load is deterministic and limit states were pushed as far away from the 
baseline design as possible in order to demonstrate the accuracy and limitations of code 
RESTRUCT. 

64-element 216-d.o.f. nlate : The clamped-ffee-free-ffee plate in Figure 6 is subjected to a uniform 
load p and has the dimensions and properties shown in the figure. Its response is computed by 
MSC/NASTRAN and redesign is performed by RESTRUCT. The incremental optimization 
problem is nonlinear and solved by sequential quadratic programming [7] because [Ak] is a cubic 
expression of the a e s which represent fractional changes of the plate thickness [17, 22]. The 
plate is subdivided into 8 structural groups each containing 8 finite elements. Results of redesign 
are summarized in Table 5 and show very high accuracy even for changes by a factor of 2 in 
eigenvalues and maximum deflection. 

144-element 896-d.o.f cylindrical shell : The simply supported shell shown in Figure 7 is 
subjected to hydrostatic pressure load p due to 286 meters submergence in salt water [23]. 
Dimensions [25] and properties are also shown in the figure. Its modal dynamic and static 
deflection response is computed by MSC/NASTRAN [16, 18]. The optimization problem in each 
increment is nonlinear and solved by sequential quadratic programming [7], The cylindrical shell 
I is subdivided into 5 structural groups and even though symmetry is not forced by linking ) 
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symmetric groups (1 and 5, 2 and 4) as was done in the plate redesign problem, symmetry was ] 
preserved in the redesign process. Results of code RESTRUCT are summarized in Table 6 and 
show good accuracy even for changes by a factor of 2 in eigenvalues and deflection. 

In all of the above applications, the LEAP algorithm in RESTRUCT can be pushed further by 
taking additional incremental steps if higher errors are considered acceptable. For higher accuracy, 
however, one more FE analysis may be used after about 10 increments. 

CONCLUDING REMARKS 

Several two-state problems in structural analysis, design, and redesign can be formulated by 
PAR (Perturbation Approach to Redesign) and solved by a LEAP (LargE Admissible Perturbation) 
algorithm. The objective structural design is found incrementally without trial and error or 
repeated FEAs for differences in response from the baseline design of the order of 100% or more. 
In structural reliability, PAR provides an attractive alternative to Stochastic Finite Elements and the 
Systems approach. 

Computer code RESTRUCT which implements the large admissible perturbation 
methodology, is being developed since 1983, has been tested thoroughly and has generated 
confidence in its potential to solve two-state problems. Several theoretical and numerical 
developments are under way. New types of finite elements are being introduced; new structures 
are being redesigned, such as stiffened plates and shells; new two-state problems are studied, e.g. 
submarine acoustic noise reduction, redesign for buckling objectives, redesign for stress 
objectives; a perturbation approach to reliability analysis and design is being developed; larger 
scale structures are being redesigned by postprocessing FEA results by MSC/NASTRAN.V66 
which has superelement capability. For that purpose, a supercomputer version of RESTRUCT 
running on the San Diego supercomputer has been developed. 
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by PAR (rectangular blocks Indicate problems that 
can be solved by LEAP theory , shaded blocks Indicate 
problems already solved in some form) 



Figure 3. First and second order methods for 
reliability computations 



Figure 4. 10-element. 48-d.o.f. clomped hinged beom 
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1 


Figure 5. 104-element, 192-d.o.f offshore tower 



Properties : 

E = 2.07*1 0 5 MPa 
p - 7.833*1 0*Nsec 2 /mm 4 
v = 0.3 

Length = 15.088 m 
Radius r 5.029 m 
Thickness = 5.08 cm 
Hydrostatic pressure p= 2.967 MPa 


Response : 
f| = T 2 = 13.0 Hz 
f 3 = f 4 = M 3 Hz 
f 5 = f 6 = 17.5 HZ 


Figure 7. 144-element. 896-d.o.f. cylindrical shell 


Table 1. Redesign and reliability of clamped-hinged beam 


Cam 


v 7 / v 7 

w; / w 7 

F.S.* ReAnalyna Error(%) 

0 

CPU 

(msec) 

n r 


ZJ 

# 

1 

2 

3 

4 

5 

6 
7 

r .S. HeAiiAiyui fc,rror(7<>j 

0.4588 0.4578 -0.205 

0.4588 0.4570 -0.401 

0.4588 0.4554 -0.736 

0.4588 0.4545 -0.931 

2.0358 2.0440 0.405 

2.0358 2.1124 3.762 

2.0358 2.1125 3.769 

2.0358 2.1893 7.542 

2.0358 2.0458 0.489 

2.0358 2.0568 1.033 

2.0358 2.0701 1.664 

2.0358 2.0737 1.863 

7.65 

11.40 

7.39 

20.83 

7.36 

15.09 

18.03 

31050 

52620 

48378 

75520 

68486 

90435 

105928 

9 
10 

9 

9 

11 

9 

10 

8~~ 

12 

8 

12 

10 

9 

12 

21 

11 

12 

13 

14 

15 

16 
17 

0.3 0.2892 -3.586 

0.3 0.2959 -1.379 

0.3 0.2841 -5.271 

0.3 0.2940 -2.002 

3.0 3.0122 0.407 

3.0 3.1554 5.179 

3.0 3.2042 6.801 

3.0 3.3282 10.942 

3.0 3.2152 7.173 

3.0 3.0997 3.325 

3.0 3.2969 9.898 

3.0 3.1003 3.342 

13.02 

19.52 

1131 

39.65 

15.63 

25.40 

33.59 

81131 

80156 

108912 

116314 

99915 

101394 

167822 

11 

10 

11 

9 
11 

7 

10 

10 

12 

12 

12 

10 

5 

12 

21 


4. 7, 13, 15, and 17 are aolved by the generalized invene algorithm. *F5. _ Failure State 
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Table 2. 10-element , 48-d.o.f. clamped-hinged beam : Table 3. 104-elemenl, 192-d.o.f. offshore tower structural 

' structural groups and redesign variables groups, redesign variables, and dimensions 



Table 4. Redesign and reliability of offshore tower 


Case 

# 

F.S.* 

/ l i 2 

Wi /W ! 

Re analysis 

Error(%) 

F.S.* 

TT, 2 

U»3 /W 3 

Reanalysis 

Error(%) 

0 

CPU 

(msec) 

n r 

Ila 

1 

0.6598 

0.6531 

-1.018 

— 

— 

— 

4.43 

973814 

mm 

8 

2 

— 

— 

— 

0.6598 

0.6530 

-1.030 

3.8 

925711 

H 

8 

3 

0.6598 

0.6547 

-0.786 

0.6598 

0.6541 

-0.871 

9.37 

985832 

18 

8 

11 

0.5000 

0.4871 

-2.572 

— 

— 

_ 

8.03 

1589439 

El 

8 

12 

— 

— 

— 

0.5000 

0.4844 

-3.112 

21.47 

1529708 

H 

8 

13 

0.5000 

0.4895 

-2.100 

0.5000 

0.4877 

-2.462 

13.7 

1617425 

18 

8 


p = 12; cases 12 and 13 are solved by the generalized inverse algorithm. *F.S.= Failure State 


Table 5. Redesign of 64-element 216-d.o.f. plate 



-4.479 


In all cases, n r = 7,n 0 = 5 , p = 8. 


Table 6. Redesign of Simply Supported Cylindrical Shell 


























































